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We investigate a potential of discovering lepton flavor violation (LFV) at the Large Hadron 
Collider. A sizeable LFV in low energy supersymmetry can be induced by massive right-handed 
neutrinos, which can explain neutrino oscillations via the seesaw mechanism. We investigate a 
scenario where the distribution of an invariant mass of two hadronically decaying taus (t^t^) from 
X2 decays is the same in events with or without LFV. We first develop a transfer function using this 
ditau mass distribution to model the shape of the non-LFV t^h invariant mass. We then show the 
feasibility of extracting the LFV r^/i signal. The proposed technique can also be applied for a LFV 
The search. 
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I. INTRODUCTION 

Supersymmetry (SUSY) [1 is one of the most promis- 
ing candidates of physics beyond the standard model 
(SM). In SUSY, the gauge couplings can unify at a high 
scale, which leads to a successful realization of grand uni- 
fied theories. Also, it can solve the gauge hierarchy prob- 
lem and yield a dark matter candidate to explain the 23% 
of the energy density of the Universe. However, the flavor 
sector of the minimal supersymmetric standard model 
can cause trouble due to the fact that SUSY breaking 
terms can induce large flavor changing neutral currents. 
The experimental constraints on flavor changing neutral 
currents require flavor degeneracy of the SUSY particles, 
especially for the first and second generations, if SUSY 
particles are lighter than around 2-3 TeV [2]. 

In order to solve the flavor problem, the simplest as- 
sumption is that the squark and slepton masses are uni- 
fied, flavor diagonal, and degenerate in the mass basis of 
quarks and leptons. The minimal version of this picture 
is the minimal supergravity framework (mSUGRA) [3], 
which is described by 4 parameters and a sign: The uni- 
versal scalar mass, mo, the universal gaugino mass, m 1 / 2 , 
the universal triliner coupling, Aq, the ratio of the vac- 
uum expectation values of the two Higgs fields, tan/3, 
and the sign of the Higgs mixing term, fi, in the super- 
potential. These parameters are specified at the grand 
unified scale Mq- In this model, however, the renormal- 
ization group equation splits the third generation from 
the other two because of the large top Yukawa coupling. 
The Cabbibo-Kobayashi-Masukawa matrix then allows 
this model to have observable flavor violation involving 
the third generation, e.g., the branching ratio of b — > s"/. 

In the leptonic sector, the neutrino oscillation data sug- 
gest that we also have a neutrino mixing matrix, the 
Maki-Nakagawa-Sakata-Pontecorvo matrix [1] . The light 
neutrino masses are usually generated via the seesaw 
mechanism, involving heavy Majorana masses for right- 
handed neutrinos, which are introduced as additions to 



the SM quarks and leptons. The precise seesaw formula 
for the light neutrino mass matrix with three generations 
is given by [5] 



M u =Ml(M R )- 1 M Dl 



(1) 



where Md is the Dirac neutrino mass matrix and Mr 
is the Majorana matrix, which consists of three right- 
handed neutrinos (u c ) that have masses at the scale vr-l 
corresponding to a new (local) B — L symmetry. The 
neutrino mixing angles in such schemes would arise as 
a joint effect from two sources: (i) mixings among the 
right-handed neutrinos present in Mr and (ii) mixings 
among different generations present in the Dirac mass 
matrix Md- The (physical) neutrino oscillation angles 
will also receive contributions from mixings among the 
charged leptons. 

The neutrino flavor mixings induced by the seesaw 
mechanism (as needed by oscillation data) can generate 
lepton flavor violation (LFV) effects. Within the SM, 
extended minimally to accommodate the seesaw mecha- 
nism, such effects are extremely small in any process due 
to a power suppression factor (1/vb~l) 2 , as required by 
the decoupling theorem. However, this situation is quite 
different if there is low energy SUSY. The suppression 
factor for LFV effects then becomes much weaker, tak- 
ing the value (1/Msusy) 2 - This can lead to observable 
LFV effects at low energies, as noted in a number of pa- 
pers [5J. The main difference is that, with low energy 
SUSY, LFV can be induced in the slepton sector, which 
can then be transferred (through one loop diagrams in- 
volving the exchange of gauginos) to the leptons sup- 
pressed only by a factor (1/Msusy) 2 - The experimental 
evidence for neutrino oscillation is thus a strong indicator 
that there might very well be LFV, assuming the validity 
of low energy SUSY. Searches for LFV processes such as 
t — > /17 and/or /i — > e"f can therefore be an important 
source of information on the v c mixings in Mr and/or 
family mixings in Md- 

Motivated by this natural occurrence of LFV due to 



neutrino oscillation, in this paper, we study the prospect 
of discovering LFV in the slepton sector at the Large 
Hadron Collider (LHC) using mass reconstruction tech- 
niques. We do not use any particular scenario to generate 
the LFV in the leptonic sector. Instead, we use a bottom 
up approach so that this study can be applied to any 
model of LFV. 

At the LHC, the squarks and gluinos should be pro- 
duced from the proton-proton collisions in any SUSY 
model if the masses of these colored objects are < 3 TeV. 
These colored objects decay via a cascade process ulti- 
mately into the lightest SUSY particle, which is the light- 
est neutralino (xi) in this paper. We consider the well 
motivated scenario where R-parity is conserved, so that 
the Xi is stable and can be the dark matter candidate. 
At the detector, the Xi escapes detection and forms miss- 
ing energy. In order to be able to see the LFV from the 
cascade decays, the cascade decay chain must contain 
sleptons. This happens in a large class of SUSY mod- 
els where the squarks and gluinos decay into %° or the 
lightest chargino, Xi , by emitting quarks (jets). The x\ 
and the Xi then decay into sleptons and leptons. The 
sleptons eventually decay into leptons and the xi- The 
entire decay process is: 



,-o 



g^qq^ qqxl -► qqU -> qqlixi 



(2) 



It is possible that Xi rnay be produced instead of x®i m 
which case Xi decays into ivx\ via charged slepton or 
sneutrino production. 

In the mSUGRA framework, the first two generation 
sleptons, e and p,, are highly degenerate without LFV 
effects. However, introducing some LFV effects causes 
a mass splitting between these SUSY particles. If these 
particles are accessible at the LHC, then measuring this 
mass splitting can help to probe the LFV [7|- How- 
ever, the first two generation sleptons are very difficult to 
probe by direct production at the LHC. Additionally, for 
larger tan/3, the staus are much lighter compared to the 
first two generation sleptons. Thus, the branching ratio 
into these heavier sleptons from other SUSY particles is 
very small. We propose a new technique in this paper 
that only depends on a mass measurement of the lightest 
slepton. 

Our final state of interest contains at least two hadron- 
ically decaying taus, t\, plus missing transverse energy, 
pr- Using the 2j + 2r h + fir signal, it has been shown in 
various studies [5l412| that the ThTh invariant mass distri- 
bution, m TT , involving opposite sign (OS)— like sign (LS) 
combinations of the two taus shows a clear endpoint or 
peak, which is a function of the lighter stau mass and 
the two lightest neutralino masses. In particular, it is 
also shown that using kinematic observables (such as var- 
ious invariant mass distributions of the jets and taus as 
well as the pr distribution of the taus), one can solve for 
the gluino, squark, lighter stau, and lightest neutralino 
masses in various SUSY models [9j . Since the masses are 
reconstructed from the observables, this technique is very 
general and can be applied to any generic SUSY model. 



Since the slepton-neutralino and slepton-chargino in- 
teractions carry the information of LFV, we investigate 
one of the major cascade decay chains involving the slep- 
tons. As we have discussed above, in non-LFV scenar- 
ios, the final states contain only taus when staus are the 
lighter of sleptons, and the m TT distribution shows a clear 
endpoint and peak. However, when LFV is present, we 
can have the decay modes x® ~~ ^ r / i X? and/or x\ ~^ Te Xi- 
The authors of [12] proposed a search for the LFV sig- 
nal in an excess in OS /u-Th over OS e-r^ events. Their 
analysis assumed no LFV decays in the e-r^ channels. In 
this paper, we propose a complementary method using a 
"transfer" function, which can be used for the [i-Th and 
e-Th LFV channels simultaneously, (see Sec. IV). 

The paper is organized as follows. In Sec. |II| we dis- 
cuss further the origins of LFV terms in SUSY, as well 
as our particular implementation of such LFV terms for 
this study. In Sec. |III| we show "measurements" of the 
SUSY masses using non-LFV decays. Measuring these 
SUSY masses is a crucial step in estimating the non- 
LFV background. In Sec. |IV] we propose a technique to 
estimate this background and extract the effects of LFV. 
This technique is based upon the hadronic tau pair signal 
that is common to both the LFV and non-LFV case. We 
report the feasibility of detecting the LFV component at 
the LHC. Finally, we present our conclusions in SeclVl 



II. ORIGIN AND IMPLEMENTATION OF LFV 

LFV effects can be generated by the neutrino seesaw 
mechanism in SUSY as follows. At the grand unified 
scale we have the mSUGRA boundary condition and 
thus, there is no flavor violation anywhere except in 
the Yukawa couplings. It should be noted that in the 
absence of neutrino masses, there is only one leptonic 
Yukawa matrix, Yj, for the charged leptons, which can 
be diagonalized at Mq and which will remain diagonal 
to the weak scale. Thus, it would not induce any fla- 
vor violation in the slepton sector in mSUGRA models. 
However, to satisfy the neutrino mixing data, the right- 
handed neutrinos have masses of order vb-l, which is 
lower than Mq. These right-handed neutrino masses 
must be ~ I0 12 — 10 15 GeV. In the momentum regime 
Vb-l < A* — Mq where the v c fields are active, the 
soft masses of the sleptons will feel the effects of LFV in 
the neutrino Yukawa sector through the renormalization 
group evolution. At the scale vb-l, the slepton mass 
matrix is no longer universal in flavor, and this nonuni- 
versality will remain down to the weak scale. 

Since we are interested in the phenomcnological as- 
pects of such LFV at the LHC, we simply introduce a 
term which causes LFV by hand within the charged slep- 
ton mass matrix. The charged slepton mass matrix is a 
6x6 matrix and is given by 



Mi = 



M\ L M\ R 



M?. B M 2 



l LR 



RR 



(3) 



where M\ L represents the 3x3 matrix for the soft masses 
for left sleptons, M. 2 RR represents the 3x3 matrix for the 
soft masses for right sleptons, and M\ R represents the 
3x3 diagonal matrix with elements m; (Ai + fi tan 0) , 
where Ai is the trilinear soft mass term and mi is the 
diagonal charged lcpton mass (YiVd) for generation I. In 
mSUGRA, 



Mir, =Ml n = m, 



RR 



1 

1 
1 



(4) 



and A[ = A . The effects of LFV can be produced when 
off-diagonal elements in M\ L LR RR are present. During 
the diagonalization of this mass matrix (which puts us 
in the mass eigenstate basis of sleptons), if there is no 
such off-diagonal LFV element in M. 2 LL LR RR , then the 
mass eigenstates are states of pure flavors. However, the 
off-diagonal LFV element causes the mass eigenstates to 
become mixtures of different flavors. These mixed-flavor 
mass eigenstates naturally act sometimes as one flavor, 
and sometimes as another. 

Consequently, if the second lightest neutralino, x%i is 
produced at the LHC from the cascade decays of squarks, 
then it can decay to tt, /u/z, or ee final states, plus fk 
due to the undetectable lightest neutralino, Xi- On the 
other hand, when LFV producing off-diagonal elements 
are introduced, the final states from the x.2 decay can 
also have t/j,, fie, or er final states plus pr- If, for ex- 
ample, the (2, 3) element [which is the same as the (3, 2) 
element] of any or all of the M\ L LR RR are nonzero then 
t\x final states plus fir will appear in the LFV x° decay. 
In this paper, we study LFV by introducing a nonzero 
value of the (2, 3) element of the M RR matrix. This new 
element also will allow the r — > ^7 decay and therefore 
is constrained, i.e., B(r -> /ry) < 4.4 x 10~ 8 [13] . We will 
define 



'RR.LFV 



[M 2 



RR1 23 



[■M RR ] 33 



(5) 



to quantify the amount of LFV. This quantity will enter 
into the LFV decay modes of neutralinos and sleptons at 
the LHC and r — > fj.'y amplitude. 



III. DETERMINING THE MASSES OF fi, x°, 

ANDx° 

As mentioned in the Introduction, we see that the LFV 
and non-LFV decay channels both involve the f\, X\i 
and X2- in fact, for LFV which is not too large to intro- 
duce appreciable change in the stau mass, it is possible 
to reconstruct the 17/i invariant mass distribution, m rM , 
whose endpoint will coincide with the end point of the 
m Tr distribution. However, the m T ^ distribution can also 
be present even when there is non-LFV, due to the lep- 
tonic tau decay, r — ► i/Dfi. This is a major background 
to the LFV signal. 



Thus, in order to understand the LFV signal at the 
LHC, we must first be able to estimate this background. 
In order to do this, we must measure the masses of the 
SUSY particles involved in both the LFV and non-LFV 
signals as accurately as possible. Thus, in this section we 
demonstrate the technique used to determine the SUSY 
particle masses involved in the LFV decay chain, which 
we will discuss in Sec. |IV| Specifically, we will determine 
the fi, Xi, an d X2 masses. 

Since the subsystem xlr^i'Xi involves the decay chain 
that is essential for the study of LFV at the LHC, we 
need to produce this subsystem in the cascade decays of 
q, g that can occur as follows: 



<1L 



ixl 



-t qr T f^ — > 



q^T+xi 



(6) 



The signal of this decay chain at the LHC is character- 
ized by high energy jets (from the squark decays), a pair 
of oppositely charged tau leptons, and a large missing en- 
ergy signal (from the lightest neutralino, which escapes 
detection) . We need to determine the X2 1 an d 'Tl masses 
from this chain. 

In order to determine the masses, we choose a model 
in mSUGRA [3 , generate the mass spectrum for that 
model using SPheno [15], simulate LHC collision events 
at <fs = 14 TcV using PYTHIA [16], and model the de- 
tector response with PGS4 [17 . Although the model 
point is excluded by the LHC experiments, we choose 
m = 250 GeV, m 1/2 = 350 GeV, A Q = 0, tan/3 = 40, 
and li > as a benchmark point for a comparison with 
our previous study for the coannihilation case [S]. The 
relevant masses x§ 1 anc ^ ^1 a ^ ^ n ^ s benchmark point are 
shown in Table [T] We stress that the technique presented 
here can be used for any SUSY model point with heav- 
ier SUSY masses as long as the X2"^i"Xi subsystem is 
present in the cascade decay chains of SUSY particles. 

As stated above, our benchmark point has already 
been ruled out at the LHC, since g ~ qx ~ 800 GeV |18j . 
However, since the q and g masses set the cross section 
of SUSY production in the analysis, we can go to any 
other point of the parameter space with larger q and g 
masses. These larger masses do not change the presence 
of the xlr^i"Xi subsystem of interest. Thus, by scaling 
the luminosity by the cross section times the branching 
ratio of q and g into xlL we can achieve the same de- 
gree of accuracies for determining the X2 1 an d f\ masses 
or the amount of LFV. For this parameter space point 
£>(<?£ -> X2I) — 0-35 and B(q R -> X2Q) — 0- However, 
we can change these branching ratios by changing the 
wino-bino content of the lightest two neutralinos. For 
instance, by departing from the gaugino mass unification 
scheme in mSUGRA, we can find a model where the X2 i s 
almost entirely bino-like. In this case, B(q R — > X2I) — 1> 
increasing the overall branching ratio of squarks into our 
decay chain of interest by a factor of about three. 

Thus, the luminosity requirement in the analysis of the 
q L chain is dictated by a(q L ,g)B(q L ,g -> x^^-j+pr-^rr, 
where <t(<Zl, g) is the production cross section of q~L an d g 
at the LHC, Aj + p T is the acceptance for the jets plus ft? 



3000 



TABLE I: Relevant mass spectrum for our chosen mSUGRA 
benchmark point: mo — 250 GeV, m.1/2 = 350 GeV, Aq = 0, 
tan/3 = 40, and fi > 0. All masses are in GeV. 



Particle True mass 






186.7 
141.5 
265.8 



system, and A TT is acceptance for the tt system. Now, 
as the cjl and g masses increase, a(qL,g) will go down. 
However, Aj + f^ and A TT can be maintained the same 
(by adjusting the cuts). Thus, to achieve the same result 
for a different a x B, the analysis technique remains the 
same; only a different amount of luminosity is needed. 

Our benchmark point described above (and shown in 
Table pi) has an overall production cross section of 6.6 pb 
according to our PYTHIA simulation. If we choose a model 
point that has a similar decay chain scenario and that has 
not already been ruled out, we can estimate how much 
more luminosity would be required as compared to our 
benchmark. For example, another mSUGRA point (with 
too = 410 GeV, to 1/2 = 750 GeV, A = 0, tan/3 = 40, 
and n > 0) has larger squark and gluino masses, yet 
still has the tt system of interest. This point has a 
cross section of 0.094 pb. Thus, we would naively ex- 
pect to require a factor of 70 more luminosity than we 
present here. As suggested above, the acceptances may 
also change. However, we can alter the cuts to maintain 
the luminosity requirement. We have shown a similar 
luminosity scaling behavior in a previous analysis |12j 

Our analysis proceeds as follows. In order to select our 
SUSY events from the background of other SUSY events 
and SM background events, we employ similar cuts as 
were used in [!7J. For us to select the event for analysis, 
it must have the following: 

• at least two hadronically decaying tau leptons with 
p^ T > 15 GeV [IH], 

• at least two jets, where the leading two jets have 
PTJetl.2 > 100 GeV, 

• missing transverse energy, fjj > 200 GeV, and 



scalar sum, /iq 



-PTJetl +PT,jet2 > 600 GeV. 



Once events have been selected in this way, we select all 
pairs of tau leptons from each event. Each pair is charac- 
terized as either LS or OS based upon the reconstructed 
charge of the taus in the pair. To remove the combina- 
toric background of incorrect combinations of taus, we 
can perform the OS— LS subtraction for any kinematical 
distribution in which we are interested. Doing this, we 
construct the following kinematical distributions: 

• the ThTh invariant mass distribution, m TT , 




100 150 200 
m Tt (GeV) 



300 



FIG. 1: The ThTh invariant mass distribution for our bench- 
mark point, shown in Table IT] A linear fit finds the endpoint 
of the distribution. This distribution represents an integrated 
luminosity of 1000 fb _1 . However, we also report the situation 
for a lower luminosity of 300 fb _1 in this paper. 



• the distribution of the transverse momentum of the 
higher p T tau, Px,r(high) ' 

• the distribution of the transverse momentum of the 
lower p T tau, p^rflow)' and 

• the distribution of the transverse momentum sum 
of the two taus, slope(p™ + ) = PT,r(high) +^T%(low) ■ 

In order to determine the fi, Xi> an d X2 masses, we 
need three independent observables. Here, we over con- 
strain the system with four observables, which allows us 
to reduce the uncertainty in the measurement. The four 
observables we choose are as follows. 

The ThTh invariant mass distribution, m TT , has a max- 
imum value for taus coming from the decay chain shown 
in Eq. Q. (We select taus from this decay chain us- 
ing the OS— LS technique described above). This maxi- 
mum value, m™ d , depends on all three masses: m™ d = 
fiimf-L, T7iy0, m.yo). A sample ThTh invariant mass distri- 
bution is shown in Fig. [TJ 

The slope of the log scale plotted pt distributions 
can also be made into observables. For instance, the 
slope of the higher px tau is a function of all three 
masses: slope (p™ T(h . gh) ) = / 2 (TO fl ,TO^o,TO^o). Also, 
the average of the slopes of the high and low taus, 
(sloped™)) = i(slope(p™ r(low) ) + slope(p™ (high) )), 

is a function of all three masses as well: (slope(px IS T )) — 
/^(to^ , m^o, m$o). Lastly, the slope of the transverse 
momentum sum distribution is a function of two of the 
masses: slope(px ls + ) — j \(m^o , to^o) . 

In order to find these functional forms, /1, f%, fa, 
and /4, we vary the masses near our benchmark point. 
We change each of the masses, rrif 1 , m^a, and to-^o, in 
turn while holding the others constant. By repeating 
the simulation for each of these varied points, we can 
find the functional forms of each observable as a func- 
tion of one mass. The overall functional forms arc then 



TABLE II: Mass measurements for our chosen mSUGRA 
benchmark point: mo = 250 GeV, mi/s = 350 GeV, Aq = 0, 
tan/3 = 40, and /j, > 0. The statistical uncertainties shown 
first are estimated for a luminosity of C = 1000(300) fb _1 . 
The systematic uncertainty due to a jet energy scale mismea- 
surement of 3% |20[ has also been estimated, shown second. 
All masses are in GeV. 



TABLE III: The effects of Srr^fv upon our benchmark point 
(shown in Table III. The LFV decay \2 ~ > A^i has a negligible 
decay branching ratio (B < 10~ 4 ), which we do not show here. 
However, this decay is still included in our simulations. We 
note that values of 5hr,lfv larger than ~ 15% violate the 
bound on the branching ratio B(j — > fi'y) < 4.4 x 10 -8 [2] 
for our benchmark model. 



Particle mass 


Solution one 


Solution two 


fi : 186.7 
X? : 141.5 
X° ■ 265.8 


181.5 ±3.7(5.1) ±4.1 

140.6 ±5.4(6.5) ±6.2 
265.3 ±6.2(8.5) ±7.3 


205.8 ±5.9(6.1) ±5.7 
151.4 ±6.4(8.6) ±6.3 
278.9 ±9.2(11.7) ±9.0 



estimated by combining these one-dimensional functions 
into a three-dimensional function in a multiplicative way. 
Similar mass determination techniques (using kinemati- 
cal observables and functional forms) have been demon- 
strated before. See for example [B1-IT2]. 

Once the functional forms are found in this way, we 
can invert them to solve for all three masses. In prin- 
ciple, this can be done algebraically. However, such a 
system of equations is complicated and over constrained. 
Instead, we invert these equations by use of the Nelder 
Mead method. This method is a nonlinear optimiza- 
tion technique that we employ to search for the masses 
that best fit the observables according to the functional 
forms. This is nearly identical to the method used in 
[12]. We find two solutions for the SUSY masses using 
this method, due to the nonlinear nature of the functional 
forms. The results of our mass determination are shown 
in Table HD 



IV. SEARCHING FOR THE LFV SIGNAL 

Now that the masses have been determined (in spite 
of having two mass solutions), we can investigate the 
effects of including the <5r.r,lfv term into our model. To 
see the effects of this at the LHC, we choose a value 
for (Sr^lfv an d use it to rediagonalize the slepton mass 
matrix for our benchmark point. The lightest slepton, 
i\ , then becomes a linear combination of \xb. in addition 
to the original tr and fx states. This allows for the LFV 
decays 



X°2 



and 



id! 



«fr 



(7) 



(8) 



Thus, the final state may include one or more muons 
from the LFV decays instead of the taus in the decay 
chain shown in Eq. Q. Thus, if we can select these LFV 
muons, we may be able to see the effect of (Jr^lfv- 

However, this analysis is complicated greatly by the 
fact that some of the taus decay naturally to muons. In 



<$RR,LFV 


(%) 


m h (GeV) B(h -> nx°i) 







186.7 





2 




186.3 


4.9 x 10~ 4 


5 




186.0 


3.1 x 10~ 3 


10 




185.1 


1.2 x 10~ 2 


15 




183.5 


2.6 x 10~ 2 



order to see the effects of <Sr.r,lfVi we need to discrim- 
inate between the muons from r decays and the muons 
from the LFV decays shown in Eqs. (l7| and ffity. 

Once we have found such a method of discrimination, 
we should be able to see the LFV signal in a kinematic 
observable similar to m TT . We plan to see the signal us- 
ing the Th/i invariant mass distribution, m Tfi . We start by 
generating the models for a few values of iJr^lfv- The 
consequences for various values of (J&r^fv are shown in 
Table |TO) In this table, we show how the lightest slepton 
mass changes (which we calculate by rediagonalizing the 
slepton mass matrix) , as well as the change in the branch- 
ing ratio for the decay shown in Eq. ph. (We calculate 
the tree-level decay width to determine the branching 
ratios.) In practice, we change the model only by intro- 
ducing the additional decay channels shown in Eqs. (fTj) 
and ^ (even though the former has a negligible branch- 
ing ratio). However, we do not bother to change the 
stau mass due to the effect of <$r.r,lfVj since it shifts 
only slightly (easily within la of our measured value in 
Table [Il|. 

We simulate the LHC signals of the models for the 
four nonzero values of g^,r,lfv using PYTHIA and PGS4 
as before. We include the branching ratios for these LFV 
decays in the input files given to PYTHIA. These four sim- 
ulations are treated as realities that we may see at the 
LHC. Thus, we treat them as "LHC data", and will refer 
to them as such. The mass determination techniques of 
Sec. Mil will result in the masses found in Table UU for 
each of these four LHC data points. 

As stated above, to find the LFV signal, we must esti- 
mate the background from muons that naturally arise 
from leptonic tau decays. To do this, we simulate a 
point with <5r.r lfv = 0. However, this simulation is 



based upon the mass determination from Sec. Ill Thus, 
instead of simulating the true benchmark point for the 
<5_r_r,lfv = case, we instead choose a point that matches 
the values for the measured masses shown in Table [TTJ In 
order to see the effect of the uncertainties shown for the 
masses in that table, we also simulate points that are 
la away from the central values of the measured masses. 
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FIG. 2: The ThTh invariant mass distribution for our LHC 
simulated 8rr,lfv = point. To generate this plot, we use 
the first solution given in Tablelll] The shape of this non-LFV 
distribution is found by the fit function given in Eq. J9). This 
distribution represents an integrated luminosity of 1000 fb~ . 
However, we also report the situation for a lower luminosity 
of 300 fb _1 in this paper. 



FIG. 3: The ThU invariant mass distribution for our LHC 
simulated 5rr,l,fv — point. To generate this plot, we use 
the first solution given in Table |H] The shape of this non-LFV 
distribution is found by the fit function given in Eq. H. This 
distribution represents an integrated luminosity of 1000 fb~ . 
However, we also report the situation for a lower luminosity 
of 300 fb _1 in this paper. 



This gives us a collection of points that we refer to as 
"LHC simulated" (Jr^lfv == points. 

First we analyze the LHC simulated (JrRjLfv = 
points to understand the shape of the m Tll distribution 
for the case of non-LFV. To do this, we form both the 
m TT distribution (as we did above in Sec. Ill ) as well as 



the m TM distribution. 



To form the 



' T jl 



distribution, we 



select events that satisfy similar cuts as above: 

• at least one hadronically decaying tau lepton with 
p T is T > 15 GeV, 

• at least one muon with Pt,^ > 20 GeV, 

• at least two jets, where the leading two jets have 
PT,jetl,2 > 100 GeV, 

• missing transverse energy, fij > 200 GeV, and 



scalar sum, /it 



-PTjctl +Pl\jct2 



> 600 GeV. 



In order to understand the shape of the non-LFV m Tll 
distribution, we relate it to the shape of the m TT distri- 
bution using a transfer junction. 

With a given small branching ratio for the LFV decay, 
we take advantage of the coexistence of non-LFV and 
LFV decays to construct the transfer function in the fol- 
lowing experimental steps: (1) We observe and measure 



the 



-.non-LFV 



shape in our LHC data point. (2) Using 



the LHC simulated (Jr^lfv == simulation (where the 
m£° n_LFV shapes in LHC data and LHC simulated are 
matched) we perform empirical fits to the m"° n ~ LFV and 
distributions. (3) The transfer function is the 



/?? 



non-LFV 



ratio of these fit functions. 

The fits that form the transfer function are shown for 
the m TT and m Tfi distributions in Figs. [2| and [3| These 
fits are used to try to minimize the statistical fluctuations 



by fitting with a "smooth" function. The empirical fit 
function we use is given by 



Km) = 



Po + Pi{m - p 3 ) + p 2 {m - p 3 ) 2 ifra<p 3 



(9) 

where m is the invariant mass that we are fitting, and 
the piS are fit parameters. We find for these distri- 
butions that the fits perform best if we fix the value 
of the transition parameter p 3 . For the T n Th invariant 
mass, we choose p 3 = 75 GeV, and for the T n /x invariant 
mass, we choose p 3 — 65 GeV. The fit range we use is 
20 GeV < m < 110 GeV. 

The transfer function, which is formed by the ratio of 
these fits, is shown in Fig. |4j For the LFV points, we use 
the transfer function to transform the m"° n_LFV distri- 

shape. Then, we subtract the 



bution into a to"™- lfv 

TO non-LFV sua p e f r0 m the m^ ta distribution. Any sig- 



"T/i 



nificant excess after this subtraction makes up our LFV 
signal. 

We use this transfer function with our four nonzero 



LHC data <5_r_r.lfv points from Table III The result of 
this transfer function method is shown in Fig. [5J This 
figure shows that this method has only a little sensitivity 
even for the largest allowed value of #rr,lfv = 0.15. This 
figure also shows that in spite of having two solutions for 
the measured SUSY masses, either mass solution will give 
us a similar result for the LFV signal. We compare the 
resulting m^ FV shape of this analysis to the expected 
number of LFV events in Fig. [6j 

We determine how significant this signal is compared to 
the null case of <5rr,lfv — in the region in Fig. [5] where 
the uncertainties propagated from the transfer function 
are not too large, namely, 60 GeV < m rp < 110 GeV. 
We do this by assuming a Gaussian distribution for each 
bin based upon these uncertainties. We find that the 
LFV signal for the first (second) SUSY mass solution 
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FIG. 4: The transfer function for our two LHC simulated 
<5rr,lfv = points, each based upon one of our two SUSY 
mass solutions shown in Table [III The first solution is the 
solid black data points, and the second solution is the dashed 
black data points (which are shifted horizontally slightly for 
visibility). Notice that the transfer function does not vary so 
much between the two SUSY mass solutions. These trans- 
fer functions are generated using an integrated luminosity of 
1000 fb -1 . However, we also report the situation for a lower 
luminosity of 300 fb" 1 in this paper. 
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FIG. 5: The t^jj, invariant mass distribution for our LHC data 
<5rh,lfv = 0.15 point. The shape of the non-LFV distribution 
is estimated by using the transfer function as described in the 
text for the mass range 20 GeV < m TM < 110 GeV. This plot 
shows the result for both of the mass solutions shown in Ta- 
ble [Hi The second solution is dashed and shifted horizontally 
slightly for visibility. Note that the uncertainties shown for 
the resulting (black) m^ v histograms include our estimate 
of the systematic uncertainties. This result represents an in- 
tegrated luminosity of 1000 fb -1 . However, we also report the 
situation for a lower luminosity of 300 fb -1 in this paper. 



shown in Tabic [TT] has a 2.2a (1.6c) excess for a lumi- 
nosity of 1000 fb - . We also note here the estimated lu- 
minosity requirements to obtain an equivalent signal for 



different values of <5r.r.lfv> which we show in Table IV 
For a luminosity of 300 fb , this significance drops to 
1.7(7 (1.2a). We note that the change in uncertainties on 
the SUSY mass and transfer function determinations do 
not scale as yC in our analysis. This indicates if the sys- 
tematic uncertainty in our determination technique can 



FIG. 6: The black histogram shows the LFV Th/U invariant 
mass distribution for our LHC data 5r.r,lfv = 0.15 point, 
compared to the filled red histogram, which is the Monte 
Carlo truth distribution of LFV decays for the same data 
set. This plot shows the result for both our mass solutions 
shown in Table ILTl The second solution is dashed and shifted 
horizontally slightly for visibility. Note that the uncertain- 
ties shown for the resulting (black) m^ v histograms include 
our estimate of the systematic uncertainties. This result rep- 
resents an integrated luminosity of 1000 fb -1 . However, we 
also report the situation for a lower luminosity of 300 fb -1 in 
this paper. 



TABLE IV: The required luminosity to obtain an equivalent 
signal for different values of 5r.r,lfv at our benchmark point. 
We estimate these luminosity values using our basic result of 
a 2.2cr (1.6a") excess of LFV signal for our benchmark model 
with 5_r_r,lfv = 0.15. We note that values of 5rr,lfv larger 
than ~ 15% violate the bound on the branching ratio B(t — > 
/jTf) <4.4x 10" 8 QI]. 



r {%) B(li => MX?) C (fb" 1 ) 



5 


3.1 x Id' 6 


8390 


10 


1.2 x 10 -2 


2170 


15 


2.6 x 10 -2 


1000 


32 


1 x 10" 1 


260 


45 


2 x 10" 1 


130 



be improved, the discovery potential is enhanced. 

If there were no systematic uncertainty for measuring 
the Th energy scale, the significance at 300 fb -1 would 
bounce back to 2.2a (1.6a). Thus, a decrease in the sys- 
tematic energy scale uncertainty at the LHC experiments 
would greatly improve the significance of this signal. 



V. CONCLUSIONS 

In this paper, we investigated the possibility of finding 
evidence of LFV using a new technique, called a "transfer 
function," along with mass reconstruction techniques at 
the LHC. We constructed the t^t^ and t^/j, invariant mass 
distributions. We used the transfer function to convert 



the ThTh mass distribution into the non-LFV Th/i mass 
distribution, rn"° n_LFV . The subtraction of the non-LFV 
Th/x distribution from the regular t^jjl distribution left us 



LFV 



with the LFV signal distribution, m^ 

For our benchmark model (<5r.r,lfv = 0.15), shown in 
Table |lj we needed ~ 300 fb^ 1 to observe a ~ la LFV 
signal for the case where B(£\ —> /UX?) — 3%. One can 
probe (at the - 2cr level) B{I X -> nxl) - 20% with 
a luminosity of ~ 45 fb _ for the parameter space dis- 
cussed so far. However, these luminosity requirements 
are highly model dependent. If one goes to any other 
model point where \2 decays into 1\t — > ttxI, our analy- 
sis still applies. The luminosity requirement can be scaled 
by aq t gB(q,g — > jq)- O ne can reduce the requirement of 
luminosity with a reduction in the systematic uncertainty 
of the tau energy scale at the LHC experiments. 

Lastly, we emphasize here that we have developed a 
new technique that can probe a LFV effect in this com- 



plex final state at the LHC. This technique will be equally 
effective to probe LFV in other SUSY models with a sim- 
ilar decay chain as a final state. 



VI. ACKNOWLEDGEMENTS 

We would like to give thanks to Sascha Bornhauser for 
his early work on this project, to Nathan Krislock for 
introducing us to the Nelder Mead method, to Kechen 
Wang and Kuver Sinha for their help in getting it to 
converge, and to Youngdo Oh for his helpful comments. 
This work is supported in part by DOE Grant No. DE- 
FG02-95ER40917 and by the World Class University 
(WCU) project through the National Research Founda- 
tion (NRF) of Korea funded by the Ministry of Educa- 
tion, Science, and Technology (Grant No. R32-2008-000- 
20001-0). 



[1] S. P. Martin, "A Supersymmetry primer," In *Kane, 
G.L. (ed.): Perspectives on supersymmetry* 1-98. hep- 
ph/9709356 ; H. P. Nilles, Phys. Rept. 110, 1-162 (1984). 

[2] F. Gabbiani and A. Masiero, Nucl. Phys. B 322, 235 
(1989); J. S. Hagelin, S. Kelley and T. Tanaka, Nucl. 
Phys. B 415, 293 (1994); F. Gabbiani, E. Gabrielli, A. 
Masiero and L. Silvestrini, Nucl. Phys. B 477, 321 (1996). 

[3] A. H. Chamseddine, R. L. Arnowitt, P. Nath, Phys. Rev. 
Lett. 49, 970 (1982); R. Barbieri, S. Ferrara, C. A. Savoy, 
Phys. Lett. B 119, 343 (1982); L. J. Hall, J. D. Lykken, 
S. Weinberg, Phys. Rev. D 27, 2359 (1983); P. Nath, R. 
L. Arnowitt, A. H. Chamseddine, Nucl. Phys. B 227, 121 
(1983). 

[4] B. Pontecorvo, Zh. Eksp. Teor. Fiz. 33, 549 (1957), 34, 
247 (1957), and 53, 1717 (1967); Z. Maki, M. Nakagawa 
and S. Sakata, Prog. Theor. Phys. 28 870 (1962); V. N. 
Gribov and B. Pontecorvo, Phys. Lett. B 28, 493 (1969). 

[5] P. Minkowski, Phys. Lett. B 67, 421 (1977); T. 
Yanagida, KEK-79-18, in proc. of KEK workshop, eds. O 
Sawada and S. Sugamoto (Tsukuba, 1979); P. Ramond 
CALT-68- 709, talk at Sanibel Symposium (1979) hep- 



ph/9809459 ; S. L. Glashow, HUTP-79-A059, lectures at 



Cargese Summer Inst. (Cargese, 1979); M. Gell-Mann, 

P. Ramond and R. Slansky, in Supergravity, eds. P. 

van Nieuwenhuizen and D. Z. Freedman (North-Holland, 

Amsterdam, 1979); R. N. Mohapatra and G Senjanovic, 

Phys. Rev. Lett. 44, 912 (1980). 
[6] F. Borzumati and A. Masiero, Phys. Rev. Lett. 57, 961 

(1986); L. Hall, V. Kostelecky and S. Raby, Nucl. Phys. 

B 267, 415 (1986); J. Hisano, T. Moroi, K. Tobe, M. Ya- 

maguchi and T. Yanagida, Phys. Lett. B 357, 579 (1995). 
[7] A. J. Buras, L. Calibbi, P. Paradisi, J. High Energy Phys 

06, 042 (2010). 
[8] I. Hinchliffe, F. E. Paige, M. D. Shapiro, J. Soderqvist, 

and W. Yao, Phys. Rev. D 55, 5520 (1997); I. Hinchliffe 

and F. E. Paige, Phys. Rev. D 61, 095011 (2000). 
[9] R. L. Arnowitt, B. Dutta, A. Gurrola, T. Kamon, 

A. Krislock, and D. Toback, Phys. Rev. Lett. 100, 231802 

(2008) 
[10] B. Dutta, A. Gurrola, T. Kamon, A. Krislock, A. B. La- 



hanas, N. E. Mavromatos, and D. V. Nanopoulos, Phys. 

Rev. D 79, 055002 (2009) 
[11] B. Dutta, T. Kamon, A. Krislock, N. Kolev, and Y. Oh, 

Phys. Rev. D 82, 115009 (2010) 
[12] B. Dutta, T. Kamon, A. Krislock, K. Sinha, and 

K. Wang, |arXiv:1112.3"966] [hep-ph]] . 
[13] E. Carquin, J. Ellis, M. E. Gomez, S. Lola, and 

J. Rodriguez-Quintero, J. High Energy Phys 0905, 026 

(2009). 
[14] K. Hayasaka et al. [Belle Collaboration], Phys. Lett. B 

666, 16 (2008); B. Aubert et al. [BABAR Collaboration], 

Phys. Rev. Lett. 104, 021802 (2010). 
[15] W. Porod, Comput. Phys. Commun. 153, 275 (2003). 

We use SPheno version 3.1.4. 
[16] T. Sjostrand, S. Mrenna, and P. Skands, J. High Energy 

Phys. 05, 026 (2006). We use PYTHIA version 6.411 with 

TAUOLA. 
[17] PGS4 is a parameterized detector simulator. We use 

version 4 (http://www.physics.ucdavis.edu/~conway/ 

research/sof tware/pgs/pgs4-general .htm) in the 



CMS detector configuration. We assume the r identifi- 
cation efficiency with p^ s > 15 GeV is 50%, while the 
probability for a jet being misidentified as a r is 1%. The 
fr-jet tagging efficiency in PGS is ~42% for E T > 50 GeV 
and |r/| < 1.0, and degrading between 1.0 < |~/| < 1.5. 
The 6-tagging fake rate for c and light quarks/gluons is 
~ 9% and 2%, respectively. 

[18] CMS Collaboration, Phys. Rev. Lett. 107 221804 (2011); 
ATLAS Collaboration, J. High Energy Phys. 11 099 
(2011). 

[19] CMS Collaboration, "Search for Physics Beyond the 
Standard Model in Events with Opposite-sign Tau Pairs 
and Missing Energy", CMS Physics Analysis Summary 
CMS-PAS-SUS-11-007 (2011). 

[20] The systematic uncertainties at the LHC will be evalu- 
ated correctly by the ATLAS and CMS collaborations. 
In this study, we demonstrate the effects of systematic 
uncertainties with a jet energy scale mismeasurement of 
3%. This mismeasurement is applied to jets, taus, and 
the missing transverse energy. 



